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Abstract 



T3 

C I The flux-line lattice in the compound ErNi2B2C, which has a tendency to 

ferromagnetic order in the a-b plane is studied with external magnetic field 

►^ direction close to the c-axis. We show the existence of an instability where 

(^ ' the direction of flux-lines spontaneously tilts away from that of the applied 

Cn ■ field near the onset of ferromagnetic order. The enhanced fluctuations in the 

Qv ' flux lattice and the square flux lattice recently observed are explained and 

C^ ' further experiments suggested. 
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Recently, it was discovered that the flux-hne lattice of the anisotropic magnetic su- 
perconductor ErNi2B2C has very unusual properties. In small angle neutron scattering 
experiments it was observed that the vortex lattice structure (with magnetic field applied 
at c-direction) is a square lattice as opposed to hexagonal as is typically found in a high 
field, type II superconductor |T[]. Moreover, for applied field H at an angle 9a ~ 1-6" away 
from the c-axis, the total field B is rotated towards the a-b plane and the angle /S.9 between 
H and total magnetic field B increases rapidly at low temperature [0 . The FWHM of the 
rocking curve am which measures the longitudinal correlation length ^l of vortex lines along 
their length is also found to increase sharply at low temperature, with qualitatively similar 
temperature dependence as the angle A^ |jl|. 

It is also observed that ErNi2B2C has a tendency of developing weak ferromagnetic 
order at low temperatures T < 2.37^ 0,0] • Such a transition cannot occur with a uni- 
form superconductive state preserved [§]-§]. In a previous paper [^j we have suggested that a 
spontaneous vortex phase occurs in this material in the ferromagnetic state. The magnetic 
properties of the material can be described by the Ginsburg-Landau free energy functional 

+ lc,|M|2 + l(3\M\^ + ^7'l VM|' - 5.M], (1) 

where 5 = V x A, M is magnetization and ip is the superconducting order parameter. The 
magnetic component M is found to be strongly anisotropic in ErNi2B2C, with magneti- 
zation M reside essentially only along the in plane easy axis in (100) and (010) directions 
0. We have shown that the unusual in-plane magnetic response of the compound can be 
explained using the GL functional (1) [Q. In this paper we shall study the out-of-plane 
magnetic responses of ErNi2B2C using the same GL functional (1). We shall assume that 
the magnetization M lies only on a-b plane and shall consider applied magnetic field making 
small angle 6a with the c-axis. The in-plane anisotropy of the magnetic component is not 
included in the GL functional (1) but shall be considered later to understand the square 



vortex lattice structure. 

The competition between magnetism and superconductivity appears in Eq. (1) as a 
Meissner effect of tlie superconducting component towards tlie internal magnetic field pro- 
duced by the magnetic component B = AnM. For systems with superconducting transition 
temperature T^ higher than the magnetic transition temperature T^, the magnetic transi- 
tion is suppressed. However, the system may go through a second order phase transition 
to a spiral phase or a first order transition to spontaneous vortex phase at a slightly lower 
temperature Tg < Tm ^^- We shall concentrate at the temperature region T > Ts(m) in 
this paper and shall study changes in magnetic response of the system as T — >■ Ts(m)- In this 
temperature range M is small and we can neglect the M^ term in the GL functional. The 
qualitative behaviour of the system at this temperature range can be most easily understood 
by considering the London limit where if) = constant and neglecting the |VMp term in F 

i- 

Writing B = Bc + B^b, where Be and Bab are the magnetic field along c direction and on 
a-b plane, respectively, we obtain after minimizing Fwith respect to M and A, M = Bab/ a, 
and A = AgV x B', where B'^ = B^, and B'ab = (1 —4:TT/a)Bab- Putting M and A back into 
F, we obtain 



F ~ d-^r 



'' ' ^Bl + il--)Bl + XliVxBf 



(2) 



where Ag = ■mc^/87re^|'?/'p is the London penetration depth for the 'pure' superconducting 
component, a is a decreasing function of temperature and the magnetic transition (in the 
absence of superconducting component) occurs at a(Tjn) = 47r. Notice that for magnetic 
field in the a-b plane, the presence of a magnetic component reduces the overall cost in 
magnetic energy of the 'pure' superconductor by a factor (1 — An/a) and also reduces the 



London penetration depth from from Aq to A = {Jl — 47r/a;)Ao 0. Notice that Eq.(2) 
can also be mapped into the problem of anisotropic superconductors by rescaling variables. 
However, in the presence of external field H, the two problems become very different as we 
shall see later (see also ref. and |]10|)- To study the magnetic response, we first consider 



a single vortex line solution in the London limit using Eq. (2). We shall assume that the 
vortex line is located in the a-c plane and makes a small angle 9^ with the c-axis. Minimizing 
the free energy, we obtain 

Ao(V X V X fi') + ^ = n^Q5{y)5{x - ztanO^), (3) 

where Hx = sin6y,ny = and Uz = cosO^. $o is the magnetic flux quantum. The equation 
can be further simplified by using the vector identity V x V x i?' = V(V.-B') — V^-B' and 
the Maxwell equation V.-B = 0. For small angle 6^, after some algebra and transforming to 
momentum space, we obtain to order 0(^^), 

:i-el/2)%s{q. + e,qx) 



BM 



1 + \lq^ + eiX^ql 
e^(!>o6{qz + e^qx) AnXl 0^ql%5{qz + e^q^) 



1 + A2g2 a (l + AV)(l + Agg2) 



TD (^ _ 47rAg Oyqxqy^oSjqz + Oyqx) /,x 

'^'^'~ a (l + AV)(l + Agg2)- l^i 

The various terms in Eq. (4) can be understood as follows: in the absence of the magnetic 
component {a -^ oo), only Bz{q) and the first term in B,j.{q) are non-zero and represents 
the magnetic field of a vortex line tilted away from the c-axis on a-c plane with small angle 
9y. In the presence of magnetic component on a-b plane, the magnetic response becomes 
anisotropic, leading to difference in penetration depth in c and a-b directions (Aq and A, 
respectively) and distortion of vortex core where a small net magnetization in a-direction 
is induced. As a result, a small magnetic dipolar field is induced in the a-b plane which is 
represented by By and the second term of B^- The energy of single vortex line ei can be 
computed using Eq. (2) and (4). In the limit (1 — Arc /a) « 1, we obtain to order 9l, 

ei = eo(l + ai^^), (5a) 

where 

Ti ,\^ X ., An, 1 AtxX^ , , , 

a^^-{--2){ln^r^-—^, 5b 



where ^ is the superconductor coherence length and eo ~ ($Q/47rAo)/ri(Ao/0 is the vortex 
hne energy when 9^ = 0. The first term in ai conies from the induced magnetic moment 
in the vortex core and from the corresponding dipolar field. This term is positive ; the 
other correction terms are negative and represent the lowering in energy from the magnetic 
component when magnetic field is in the a-b plane |^. Similar analysis as can also be 
made when the |VMp term is included in the GL functional. We find that the qualitative 
behaviour of the vortex solution is not modified except that the divergence in A~^ as T — ;> T^ 
is removed once the |VMp term is included. In particular, the London penetration depth 
saturates at a value of order A ~ (Ao^m)^^^ as T — > T^ for transition to spiral state 0], where 
^m ~ T^/ct is the coherence length of the magnetic component. 

In the limit H ~ iJ^i where density of vortices is low and interaction between vortices 
can be neglected, we may study 9^ as function of angle of applied field to c-axis 9a using Eq. 
(5). Consider the Gibb's energy functional 

G = F- Id'r^-^ 



47r ' 

where H is applied field and the total magnetic field B is obtained by minimizing G with 
respect to B. For small angle 9a and 9^,, Gibb's energy per unit volume is 

V ~ %''^^ + ^^^^) - 17^^ 2 — )• (^) 

Minimizing G with respect to 9jj, we obtain 

9n 



^^ {l + 2a,{H,,/H)y 



and 



where Hd = 47reo/"l'o is the lower critical field when the external field is along c-axis {9a = 0). 
Notice that 9^; > 9a when ai < 0. This may occur when X/C, is small enough, i.e. when 
the superconductor is not too strongly type II and when the system is close to magnetic 
instablity point T -^ Tg. 



The above analysis can be easily extended to the intermediate density regime H ~ several 
Hci where the density of vortices is of order (27rAQ)^^ and the magnetic field B is already 
more or less uniform in the superconductor. This is the regime of experimental interests [|I|. 
The magnetic field B in this case can calculated using Eq. (3), except that the right hand 
side of the equation is replaced by n$o S/f ^iv — Yn)5{.x — X„ — ztanO^), where _R„ are the 
positions of the vortices in the vortex lattice. Following procedures similar to above (see 



also ref [jTl| for details), we obtain in the small 6a limit, 



^- ~ An^ , h!w4^ , o. n - (8a) 



where Hc2 ~ *&o/(27r^^) is the upper critical field. Neglecting Logarithmic corrections, we 
obtain after minimizing G with respect to B and 9„, B r^ H — Hd + BgOl and 

where Bg ~ (47r/Q;)(l — 271 /a)H for H >> Hd and the angle-dependent magnetization 
M(9„) is 

Notice how 6^ is enhanced by the magnetic component in this case. In particular, 6y is 
always larger than 6a for H » Hd and is diverging at low temperature T < T^, in 
contrast to the low density limit where 6^ may be small than 6a and is larger than 6a only 
when sufficiently close to magnetic instability. Physically, the magnetic instability indicates 
a magnetic transition where spontaneous magnetization along a- or b- directions appears 
in the system at low enough temperature, leading to a magnetic-field assisted spontaneous 
vortex phase [Q where the vortex lattice tilts spontaneously from the c-axis below critial 
temperature 

Tsv ^Tm 7^(1 + 2ai), 

a H 

where a' = da/dT. Direct observation of this transition using neutron scattering or magnetic 

imaging techniqiues is suggested. 
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Next we consider the case 6a = and study thermal fluctuations in positions of vortex 
lines along the c-direction in the intermediate density regime. We consider the model free 
energy Fn for N vortex lines in a sample of thickness L, defined by their trajectories [r*,(-2)] 
as they traverse through sample with external magnetic field along the z-axis, 

+ E/ yD[rA^)]dz + -Y. V{\r-{z)-f,{z)\)dz (9) 

i=l -^0 ^ i^j -^0 

where Vo{r{z)) is a random pinning potential for vortex lines and V{r) ~ eo-ft'ol^/^o) is the 

interaction potential between vortex lines. In the following we shall replace V£,{r) by the 

average pinning potential VD{rj{z)) ~ i^o{rj — RjY /2 where Rj are sites on the fiux lattice. 

We shall also expand the interaction potential around it's minimum point V{r) ~ eQ{r /dY 

where d ~ {^q/ BY'"^ is the distance between vortex lattice sites. Note that dfj(z)/dz ~ 

6v{fj) in our analysis and the first two terms in Eq. (9) come from terms proportional 

to ^^ in the Gibb's energy (7). The tilt modulus e^i is chosen such that when all 6^{fj)^s 

are equal, the energy is given by Eq. (7) et{r) can be extracted from the free energy G 

of a fiux-line lattice configuration where all but one of the fiix-lines have angle 6v{i^j) = 0. 

The field configuration in this case can be calculated from an appropriate combination of 

eqn.(4) and the field configuration leading to eqn.(7) with ^ = 0. We shall assume that 

the fiux lines form a square lattice and retain only the nearst neighbor terms in the sum 

i ^ j in computing et{r). The qualitative properties of Fn do not depend sensitively on this 

approximation. We also ignore the ordinary contribution to the tilt moduli expressing the 

effect of increased line-length as it is negligible compared to the effects we are interested in. 

We obtain 

An Arc H 

eti - eo(— + 2ai + (1 )-— ), (10a) 

a a rici 

where we have neglected the logarithmic correction terms and et(r) = et25{r — d), where 
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Notice that OyjOa = {H/Hci){eo/eti). Notice also that at low temperature T < Tm,eti is 
a decreasing function of vortex density. In particular, en ^ at the spontaneous vortex 
phase transition temperature T — > Tgy. 

The thermodynamic properties of Fn can be obtained rather easily by observing that 
the variable z can be treated as an imaginary time and the model can be mapped into a 
quantum mechanical problem of coupled harmonic oscillators on a square lattice with h ^ T 
jT^ . The Hamiltonian of our effective quantum mechanical system is in fourier space, 

r Iniik) J- ^ 

where X^(A;) = A^~^/^ Ej e'^^^-^^rj^, fi = x,y and Pfi{k) = —iTd/dX^{k) is the canonical 
momentum conjugate to X^{k). The momentum dependent mass m{k) and frequency uj{k) 
are 

m{k) = eti + Aet2lik), 

and 

uj{k) 



Ko+Af^{k) 



\ m{k) 

where 7(fc) = sin'^{kxd/2) + sin'^{kyd/2). The thermodynamic properties of the flux-line 
lattice can be obtained easily from (11). We obtain 

< X^ik,z)X,ik,z') >= 5,, J^ -.ik)\^^z'\^ 

2m{k)uj{k) 

and 

2 1 V- ^ Td ,_, 

<r >=^y2 — - — - — = (12) 

N^.2m{k)uj{k) ^{K,d^)tt2 + e,ea 

where < r^ > is the mean square displace of flux lines from their equilibrium positions, 
kqc/^ ~ average pinning energy of the vortex lines per unit length. Notice that < r^ > 
remains flnite when tti ^ as long as there is a flnite pinning strength for the flux-line 
lattice. The longitudinal correlation length measured in neutron scattering experiment is 
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^L ~ ^{k -^ 0)^^ ~ Jeti/ kq, and am ~ d/S,L, which diverges as (eti)^^/^ as the instabihty 
towards spontaneous vortex phase is approached. In particular, we obtain a scahng relation 
between a„i and 9,^/9 a-, 



2 / >^a \ /Hcl^^Kod 



While Eq. (8) and (13) are in good agreement with the existing data in Ref. (1), data over a 
range of fields and lower temperatures is required to test the theory. Notice that our theory 
for thermal fiuctuations is essentially a Gaussian theory and Eq.(8) and (13) are essentially 
mean-field results. Derviation from mean-field behaviour is expected at temperature very 
close to transition temperature T^^ when the system is in the critical regime. We have 
estimated the size of the critical regime using the specific heat obtained from our theory and 
find that the critical regime is given by 



t < 



Kod^\^ f 1 



to J Kd^AC. 

where t = |1 — T/Ts^\ and AC is the specific heat jump across the transition in mean-field 
theory. The size of the critical regime is expected to be very small at intermediate density 
of vortices rf ~ Aq. 

Finally we make a plausible suggestion on the origin of square lattice structure of the 
fiux-line lattice observed in neutron scattering experiment. Another possible origin is Fermi 
surface anisotropy which leads to anisotropic shape of vortex core. As is evident from our dis- 
cussions, the enhancement of fiuctuations in positions of vortex lines is a result of enhanced 
magnetic response arising from the magnetic component in the system. In ErNi2B2C com- 
pound, the magnetic component is strongly anisotropic in a-b plane, with magnetization M 
reside essentially only along the a- or b- axes. Thus we expect that in a treatment of the 
system with anisotropy effect included, the fiuctuations in vortex line positions, and corre- 
spondingly, the mean-square displacement of vortex lines, < ViVj > will also be anisotropic 
and enhanced more in a- and b- directions. Effectively, the shape of the vortex core becomes 
" cross-like" , with vortex core elongated along a- and b- directions. When density of vortices 



is high enough, the repulsion between vortex cores becomes important in determining the 
vortex lattice structure. In particular, because of the small energy difference between square 
and triangular lattices |11[], it is reasonable that the vortices will arrange in a square lattice 
with angle 7r/4 with respect to a- and b- axes to maximize distances between each other. 
This structure is indeed what is observed experimentally |^. 

Notice however that the anisotropy in magnetic response is absent in the linear response 
regime and comes in only through the M^ term in the GL functional. Thus the anisotropy 
effect just discussed is expected to be strong only at temperature regime T < Tm, when 
the magnetic linear response is close to diverging and the M^ term becomes important. 
Experimentally, a square vortex lattice is observed at temperature well above the instability 
temperature determined by the divergence in 9y/9a or o"m- To explain this we suggest the 
interesting possibility that the anisotropy effect are strong even at T — T^ comparable to 
or larger than Tm if the size of the vortex core, ^ is microscopically large (the size of vortex 
core ~ ^, is found to be 135A |]I| in ErNi2B2C compound). In this case, the behaviour 
of vortex core can be described by GL functional (1) with ■?/' ~ and has an magnetic 
instability at T = T^, where T^ > T^ when terms of form r^M^I'j/'p which describes effect 
of conduction-electron polarization on superconductivity are present in the system [^. In 
this case, spontaneous magnetization occurs inside vortex cores at temperature T > Tm, 



pointing randomly in the ±a and ±6 directions [0, leading to anisotropic fluctuation effect 
observed well above Tm- 

The square lattice structure is stabilized only when the vortex-vortex interaction is sen- 
sitive to the core-shape. This can happen only for vortex density above a critical value. 
Indeed, the square lattice is found only above a critical field below which it has the conven- 
tional hexagonal shape. 

This line of argument suggests also that the square fiux-line lattice will disappear and 
replaced by a regular triangular lattice if the external field is applied along (010) or (100) 
directions. It also predicts that the triangular lattice will become distorted if the external 
field is applied on (110) direction. In the later case the vortex core will be effectively 
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elongated along a-b plane. Correspondingly, the distance between vortices joint by a line 
in a-b plane will be larger than distance between vortices joint otherwise. Experimental 
observation of flux-line lattice with applied field in (100) or (110) directions are suggested. 
We wish to acknowledge very useful discussions with Peter Ganiniel,D.Huse and U. 
Yaron. T.K. Ng acknowledge partial support by UGC Hong Kong, through RGC Grant 
no.UST623/95P. 



11 



[1] 

[2] 
[3] 
[4] 
[5] 
[6] 
[7 



[9 
[10 

[11 



REFERENCES 

U. Yaron et.al, Nature, 82, 236 (1996). 

B.K. Cho et.al., Phys. Rev. B 52, 3684 (1995). 

P.C. Canfield, S.L. Bud'ko, and B.K. Cho, Physica, C262, 249 (1996). 

T.K. Ng and CM. Varma,Phys. Rev. Lett. 78, 330 (1997). 

E.I. Blount and CM. Varma, Phys. Rev. Lett. 42, 1079 (1979). 

H.S. Greenside, E.L Blount and CM. Varma, Phys. Rev. Lett. 46, 49 (1981) 

M. Tachiki et.al, Sol. State Comm.31, 927 (1979); ibid34, 19 (1980). 

CG. Kuper, M. Revzen and A. Ron, Phys. Rev. Lett. 44, 1545 (1980). 

D. Feinberg and C Villard, Phys. Rev. Lett. 65, 919 (1990). 

G. Blatter et.al, Rev. Mod. Phys. 66, 1125 (1994). 

see for example, M. Tinkham, Introduction to superconductivity, (McGraw Hill 1975). 



[12] D.R. Nelson, Phys. Rev. Lett. 60, 1973 (1988); see also D.R. Nelson, in The Vortex State, 
edited by N. Bontemps, Y. Bruynseraede, G. Deutscher and A. Kapitulnik (Kluwer 
Academic Publishers, 1994) for a review. 

[13] Notice that magnetization residing in core of a vortex-line forms an effective one- 
dimensional system which is always disordered at finite temperature T > 0. 



12 



